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Abstract. We present a theoretical study of an ensemble of X-like 4-level atoms 
placed in an optical cavity driven by a linearly polarized field. We show that the self- 
rotation (SR) process leads to polarization switching (PS). Below the PS threshold, 
both the mean field mode and the orthogonal vacuum mode are squeezed. We provide 
a simple analysis of the phenomena responsible for the squeezing and trace the origin 
of vacuum squeezing not to SR, but to crossed Kerr effect induced by the mean field. 
Last, we show that this vacuum squeezing can be interpreted as polarization squeezing. 

PACS numbers: 42.50.Lc, 42.65.Pc, 42.50.Dv 
1. Introduction 

The principal limit in high precision measurements and optics communication is given 
by the quantum fluctuations of light. For several years, in order to beat the standard 
quantum limit, a number of methods consisting in generating squeezed states of light 
have been developed [Tj. In connection with quantum information technology the 
quantum features of the polarization of light has raised a lot of attention. The generation 
of polarization squeezing has been achieved experimentally by mixing an OPO-produced 
squeezed vacuum with a coherent field j2 El > or more recently by mixing two independent 
OPA-originated squeezed beams on a polarizing beamsplitter 0j. Several schemes 
using Kerr-like media have also been proposed jSl El El, and very recently, Matsko 
et al. proposed to propagate a linearly polarized field through a self-rotative atomic 
medium to produce vacuum squeezing on the orthogonal polarization The Kerr- 
like interaction between cold cesium atoms placed in a high finesse optical cavity and a 
circularly polarized field has been studied in our group and a field noise reduction of 40% 
has been obtained [TO] - We recently observed experimental evidence of polarization 
squeezing when the incoming polarization is linear [TT] . In this paper, we present a 
theoretical investigation of polarization squeezing generated by an ensemble of X-like 
4-level atoms illuminated by a linearly polarized field. To be as realistic as possible, 
the experimental parameters values of Ref [HI EH E] are taken as references. In the 
first part of the paper, we give a detailed study of the steady state and show that self- 
rotation is responsible for polarization switching and saturation leads to tristability. We 
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derive simple analytical criteria for the existence of elliptically polarized solutions and 
the stability of the linearly polarized solution. This steady state study is essential to 
figure out the interesting working points for squeezing. In the second part, we focus on 
the case in which the polarization remains linear (below the PS threshold) and show that 
both the linearly polarized field mode and the orthogonal vacuum mode are squeezed. 
Analytical spectra are derived in the low saturation limit and enable a clear discussion of 
the physical effects responsible for polarization squeezing; in particular, we demonstrate 
that self-rotation is associated to strong atomic noise terms preventing vacuum squeezing 
at low frequency. On the other hand, saturation accounts for the squeezing on the mean 
field and crossed-Kerr effect enables to retrieve vacuum squeezing at high frequency. 
The analytical results are compared with a full quantum calculation. Finally, we derive 
the Stokes parameters and relate their fluctuations to those of the vacuum field. The 
vacuum squeezing obtained is then equivalent to the squeezing of one Stokes parameter, 
the so-called polarization squeezing [T3] . 



2. The model 




Figure 1. X-like 4-level configuration. 



The system considered in this paper is a set of N 4-level cold atoms interacting in 
an optical cavity driven by a linearly polarized field as represented in Fig^ We denote 
A± the slowly-varying envelope operators associated with the a± components of the 
light [14J. They are defined from the standard linear polarization components 

_ A x — iA y _ A x + iA y 

A+ ~ 7T~ ' A ~ ~ ~7T~ (1) 

The atomic frequencies are both equal to ui at = = U24. The field frequency is u 
and the detunings from atomic resonance are equal on both transitions to A = u at — to. 
The 4-level system is described using collective operators for the N atoms of the 
ensemble, the optical dipoles being defined in the rotating frame associated to the laser 

N 

frequency (e.g. cri 4 = e ^*|l)i(4|j). The coupling constant between the atoms and the 
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field is defined by g = £od/h, where d is the atomic dipole and £ = \/hu/2e Sc. With 
this definition, the mean square value of the field is expressed in number of photons 
per second. As in Fig HI the population of level 3 decays with rate 2^± on level 1 and 
with rate 277 on level 2, the dipole decay rate being 7 = 77 + 7j_. We consider the 
case of saturated optical pumping and neglect the relaxation rate of the ground states 
populations. This approximation is well verified for alkali cold atoms [T^]. With these 
conventions, the atom-field hamiltonian is 

H = hg[A + a 41 + A ] + a 1A + A_a 32 + Ala 23 ] (2) 

The atomic evolution is then governed by a set of quantum Heisenberg-Langevin 
equations 



- (7 + iA)a u - igA + (a n - cr 44 ) + F 14 (3) 

- (7 + iA)a 2 3 - igA-(a 2 2 - cr 33 ) + F 23 (4) 
27x^33 + 27//CT44 - ig(A\_a u - A + a 41 ) + F u (5) 
27//^33 + 27^0-44 - ig(Ala 23 - A_a 32 ) + F 22 (6) 

- 270-33 + ig(Ala 23 - A_a 32 ) + F 33 (7) 

- 270-44 + ig(Al<j u - A + o- 41 ) + F u (8) 

Note that we have not reproduced all the atomic equations, but only those of interest 
for the following. The Langevin operators F^ are 5-correlated and their correlation 
functions are calculated via the quantum regression theorem (T^j. We consider a ring 
cavity with T the transmission of the cavity coupling mirror, uj c the cavity resonance 
frequency closest to uj and r the cavity round-trip time. The cavity dephasing is 
$ c = (oj — uj c )t. The incoming quantum fields are A± and the field equations read 

(T/2 + t$ c )A + - iga u + VTA™ (9) 
(T/2 + i§ c )A_ - iga 23 + VfA m (10) 

3. Steady-state 

3.1. Atomic steady state 

The atomic steady state is readily obtained by setting the time derivatives to zero 
and using the fact that a Langevin operator mean value is zero. Defining saturation 
parameters s± for both polarizations, 



rf0"i4 

dt 
da 23 

dt 
dan 

dt 
da 22 

~df 

da 33 

dt 

<i0"44 

dt 



dA + 

~dT 
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2g*\(A ± )\* 0|/2 
s± A 2 + 7 2 A 2 + 7 2 ' 1 J 

the atomic steady state is given by 

-igN(A + ) s 1 -igN(A-) s + 1 

<^14 = —TT ; : c , (cr 23 ) = —r-r ; g (12) 

7 + «A s + + s_ 1 + S 7 + il\ s + + s_ 1 + S 

I \ AT S ~ 1 + S +/ 2 I \ AT S + 1 + S '/ 2 f-,o\ 

(a u ) = N , (or 2 2) = N 13 

s + + s_ 1 + S s + + s_ 1 + S 

(<T33> = (<T44> = JYTS : = S+ + S ^ ' 

Cl± are the Rabi frequencies and 5 is the coupling saturation parameter which plays 
a symmetrical role with respect to both polarization components. For an x-polarized 
field, S = s + = s_ = s x /2 is directly related to the intracavity field intensity 

3.2. Polarization switching 

It is well known that such a coupled system may exhibit polarization switching when 
driven by a linearly polarized field ^3 E|. in fact, the intracavity field intensities 
depend on the atomic dephasings $± and absorptions a± 

$ ± = 2$ S 2 : ~ g , a± = 2a g -±- (15) 

S + + S_ 1 + S S + + S_ 1 + D 

^ 2 A A7/7 
$ °-2(A 2 + 7 2 ) ' 00 " 2(A 2 + 7 2 ) (16) 

with $o an d ao the linear dephasing and absorption in the absence of saturation. 
These quantities depend in turn on the intensities to yield a complex coupled system. 
In order to derive analytical criteria for polarization switching, we follow the method 
given in [IS] and decompose dephasings and losses into their linear and non-linear parts, 

$± = $ z ± <5>sr with $j = — — , $ SR = Qtxsn (17) 

1 + D 

a± = ai ± a SR with a t = - , a SR = aiX S R (18) 

where &sr and asR are the non-linear circular birefringence and dichroism, related 
to the ellipticity e [T9"] 

x S r = S ~ , g+ = ~ sin 2e (19) 

s + + s_ 

Thus, as pointed out in the literature [2011211, the optical pumping induces non- 
linear self-rotation (SR) of elliptically polarized light. It will be shown in the next 
section that this effect is responsible for PS in a cavity configuration. Let us first focus 
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on the solution for the a± components. Normalizing all the dephasings and absorptions 
by T/2 (Sj = 2$j/T and aj = 2a i /T), Eqs ©,(1101) read in steady state 

S± = (l + a^asny + iS^SsR-Sc)* (20) 

with s max = 2/Ts l x the maximal intracavity intensity in the absence of absorption. 
Replacing (|20|) in (j!9|) . we derive the equation for xsr- non zero solutions correspond 
to elliptically polarized states. After straightforward calculations, we obtain 



x SR = or ( a * + 5?y SR = 6? + a*-5 2 c -l (21) 

The first trivial solution corresponds to the linearly polarized field. It follows from 
the second equation and (fT7 )) .(fT8 |) that elliptically polarized states may exist as soon as 
the existence criterion C ex is satisfied 

^ = |^-5 C 2 -1>0 (s + ^s_) (22) 

Note that the absorption brings a positive contribution to the existence of 
asymmetrical solutions: this is due to the fact that non-linear circular dichroism 
produces " self-elliptization" of the field. However, this criterion gives no information on 
the stability of the solutions. In order to get some physical insight into this complicated 
problem it is useful to look at the evolution of the linearly polarized solution. 



3.3. Interpretation of polarization switching 

In this section, we give a simple interpretation of PS as the threshold for laser 
oscillations. Let us consider the linearly polarized solution along the x axis. The 
adiabatic elimination of the atomic variables leads to 



1 dA 2 

--tt = -(1 + tic)Ay + {iSi - a t )A y - (S SR + ia SR )A x + — Af (23) 



where k = T/2r is the intracavity field decay rate. In (j2*B^) all terms have zero mean 
value and are of order 1 in fluctuations ((A x ) ^ 0). Using xsrA x = i(A y — A^/\A x \ 2 A^ y ), 
one obtains 



1 HA A 2 9 

-~df = -(! + tfcK + («« - a ^]Aj A t + ^f A 7 ( 24 ) 

Owing to SR the fluctuations of the orthogonal mode undergo a phase dependent 
gain. A similar equation has already been derived in previous theoretical works in a 
single pass scheme 0. In our configuration the presence of the cavity will lead to 
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Figure 2. Resonance curves in the optical pumping regime. The parameters are 
So = 5, s max = 0.01. The cavity dephasing corresponding to PS is Sps = 4.85, close 
to Sq as given by the PS criterion. When the cavity detuning is scanned from the 
right, the linear solution is stable until 5 C = Sps and unstable afterwards. Then the 
elliptically polarized solutions, ,s+ and s_, become stable for |i5 c | < Sps- We plot also 
the resonance peaks (dashed line) for the cavity with or N atoms, in the absence of 
SR phenomenon. 



oscillations of this mode as soon as the phase sensitive gain is larger than the losses. 
This condition may be expressed as follows 



Obviously, the linearly polarized solution is not stable when Cps > 0. However, 
the adiabatical elimination of the atomic variables does not a priori take all causes for 
instability into account. Yet, we checked that this threshold analysis was consistent 
with a numerical calculation of the atom-field stability matrix. In the following we use 
Cps as a stability criterion for the linearly polarized solution. Nevertheless, it does not 
yield information on the stability of the elliptically polarized solutions, which has been 
evaluated numerically. 

Besides, the ability of a system to produce squeezing being closely related to its static 
properties, the fluctuations of the vacuum field are expected to be strongly modified 
in the vicinity of the PS threshold. Since Eq is similar to that of a degenerate 
optical parametric oscillator (OPO) below the threshold [22], perfect squeezing could be 
obtained via SR. However, the atomic noise is not included in (I24J) and is to be carefully 
evaluated. 

3.4- Optical pumping regime 

PS is caused by a competition between the two a± optical pumping processes. We can 
understand the main features of this effect by restraining ourselves to the case where 
absorption and saturation are negligible: A 3> 7 and Sj < 1. Neglecting the excited 
state populations, the optical pumping equations for the ground state populations are 




-5 C 2 -1>0 



(25) 
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dt 

da 2 2 
~df 



= -J±s + a u + j±s_a 22 (26) 

pumping 

= ~lLS-a 2 2 + 7-L s +°"ii (27) 

pumping 



so that the a + component tends to pump the atoms into level 2, the <r_ into 1, and, 
in steady state, an oc s_ and 022 oc s + . The circular birefringence Ssr is proportional 
to the ground state population difference, and consequently, to the intensity difference 
s + — s_ [see 1)19)1] . This simple analysis allows for relating self-rotation to competitive 
optical pumping and will help us interpret the resonance curves. 

Under the previous conditions both criteria ()22|) and (|23|) are equivalent and it follows 
that the linearly polarized solution bifurcates into an elliptically polarized state for 
\$c\ < dps — \/^o ~~ 1- Consequently, PS is observed as soon as the linear dephasing is 
greater than half the cavity bandwidth (5q > 1). This represents an easily accessible 
condition from an experimental point of view : in our cesium experiment using a 
magneto-optical trap jniGHj, the number of atoms interacting with the light is iV ~ 7 10 6 . 
To find realistic experimental parameters, we assimilate each one of our X-model 
transitions to the transition 6Si/ 2 (F = 4) — 6P 3 / 2 (F = 5) of the D 2 line of 133 Cs, 
for which / y/2n = 2.6 MHz. The square of the coupling constant g is proportional to 
the ratio of the diffusion section at resonance to the transversal surface S — 0.1 mm 2 
of the beam, g 2 = ^X 2 /AttS = 4.24 Hz. The cavity transmission is 10%. To obtain a 
sufficiently high non-linearity, keeping the absorption low, a good detuning is A ~ 2O7, 
so that an approximate value for the linear detuning is 5o ~ 5. Note that the saturation 
parameters of (fTTj) are simply 

s± - At (28) 

The saturation intensity being I sat = e cy 2 h 2 /d 2 = 1.05 mW/cm 2 0, typical values 
for s± are 0.1-1. 

In Fig El are represented the admissible intensities for the a + and a_ components 
versus the cavity detuning for typical experimental values of the parameters. The peak 
centered on 5 C = 5q corresponds to the symmetrical solution. When the cavity is scanned 
from right to left, the linearly polarized field (s + = s_) intensity increases until the PS 
threshold is reached (5 C = 5ps — S ). Then one elliptically polarized state becomes 
stable. The predominant circular component, say a + , creates, via the optical pumping 
process (|26H27|) . a positive orientation of the medium an — 0, a 22 ~ N. Since the atomic 
dephasing decreases to zero for the a + component (5 + — 0), as if it were propagating 
in an empty cavity. Hence, the solution draws close to the zero-dephasing peak, that 
is, close to resonance in the range \S C \ < 5ps- On the other hand the cr_ component 
"sees" all the atoms (<5_ ~ 2<5 ) and breaks down to fit the peak centered on S c = 25 , 
which is far from resonance. In order to illustrate this interpretation of the resonance 
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Figure 3. Upper plot: resonance curves for <5o = 7, A = 2O7 and s ma x = 2. The 
linear absorption is ao = 0.35. Dashed parts indicate unstable solutions. The switching 
occurs for Sps = 2.6 and Sps = 1-5. C ex = for 5 ex = 5.6, so that the tristability 
range is 2.6 < S c < 5.6. The arrows on the hysteresis cycle correspond to increasing 
and decreasing cavity detuning scan. Below are plotted the two criteria: Cps giving 
the stability of the linear solution (plain) and C ex giving the existence of asymmetrical 
solutions (dashed). 



curves, the two Airy peaks centered on 5 C = and 5 C = 25q are represented in Fig 
121 As the cavity detuning is decreased both asymmetrical solutions reunite when the 
criterion (J25)l is no longer satisfied and the linear solution becomes stable again. These 
simple interpretations will help us understand the much more complex general case, 
when absorption and saturation come into play. 

As discussed in Ref [IB], taking into account the ground state relaxation rate 70 yields 
tristability in the unsaturated optical pumping regime (j±s± <C 70) • We will now show 
that the optical saturation also leads to tristability. 

3. 5. Tristability 

It is well-known that saturation may induce multistability for the linearly polarized 
field in our configuration and substantially modify the steady state. When the 
non-linearity is sufficient, there may be three possible values for the x-polarized field 
intensity. Therefore, saturation is an additional cause of instability for the symmetrical 
solution. In Fig El we plotted the same curves as in FigEl but for higher values of s max 
and 5q. As expected, the linearly polarized state solution is distorted as a consequence 
of the non-linear effect. The effect of absorption is also clear: whereas the symmetrical 
peak height is reduced, the a + dominant peak height is not. Indeed, the a + component 
"sees" no atoms after the switching. 

Besides, the system now exhibits tristability for a certain range of the cavity detuning. 
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Figure 4. Bistability curve: linearly polarized field intensity S = s x /2 as a function of 
Smax- PS is the switching threshold, HT and LT the higher and lower turning points. 
For certain incident intensities s ma x three solutions (1,2,3) exist for the intracavity 
intensity s x /2, of which only one (3) is stable, instead of the usual two (1,3) in the 
absence of PS phenomenon. 



As mentioned previously the existence of asymmetrical solutions is related to the pos- 
itivity of C ex , whereas the stability of the symmetrical solution is given by Cps- For 
instance, on Fig El C ex = for S c = 5 ex = 5.6 and the threshold Cps = is reached 
for a dephasing 5 C = 5ps = 2.6. Thus, in the range 5p$ < 5 C < 5 ex , two different 
sets of asymmetrical solutions exist, in addition to the linear polarization state. This 
phenomenon is due to the saturation experienced by the o~ + and cr_ components. As 
expected, we checked that only the lower branch of each asymmetrical curve is stable, 
leading to tristability for the polarization state : linear, cr + -dominant or a_-dominant. 
The system switches for a different value of the cavity detuning if the cavity is scanned 
from left to right, or from right to left (see FigEJ). Hence, unlike the unsaturated case, 
saturation induces a multistable behavior and a hysteresis cycle now appears in the 
resonance curve. 

This brief study of the resonance curve for typical parameters leads to an essential 
observation: the lower branch of the bistability curve for the linear polarized field is 
not stable. We may wonder if there is a domain of the parameter space for which it 
is not the case. Since the quantum fluctuations are expected to be most reduced in 
the vicinity of the lower turning point the answer is of crucial importance for 

squeezing and will be treated in the next section. 

3. 6. Competition between SR and saturation 

To complete the analysis of the steady state, we would like to emphasize that, when 
the cavity is scanned, PS always happens before reaching the higher turning point of 
the bistability curve. In order to get some insight into this complicated problem, it 




Figure 5. Bistability curves S = s x /2 as a function of s max for field A x . The 
three S-shaped curves correspond to different values of the cavity dephasing (S c — 
—0.25, 0.3, 1.1 from left to right), <5o and A having the same value as in FigEI The AB 
segments represent the HT (higher) and LT (lower) curves. The dashed curve is the 
ensemble of the intensities Sps for which polarization switching occurs. The system 
exhibits bistability for sa < s rnax < sb- 

is worth looking at Fig EJ We plotted the typical S-shaped variation of the linearly 
polarized field intensity s x versus the incoming intensity s max , for a fixed value of the 
cavity detuning 8 C . We choose the parameters so that there is bistability for this state 
of polarization and report the position of the lower (LT) and the higher (HT) turning 
points. In the absence of the PS phenomenon, the solutions between HT and LT are 
unstable (like 2 on FigEJ), whereas solutions on the lower (1) and higher (3) branches are 
stable. However the stability of the linear polarization is modified by the PS effects. To 
a fixed value of the dephasing corresponds the PS intensity Sps cancelling Cps in (J2HJ); 
if S = s x /2 < Sps(8 c ), then the linear polarization is unstable. Hence, if Sps(S c ) > Sht 
is satisfied in the whole parameter space, then PS occurs before reaching HT, and, 
consequently, the lower branch is never stable. 

This general feature is shown on Fig in which we represented different bistability 
curves as in Fig 0] The upper branch of AB is the HT curve (the ensemble of the 
higher turning points when S c is varied), the lower branch is the LT curve. The dashed 
curve shows the ensemble of the intensities Sps for which the polarization switches. 
This curve is always above the HT curve, confirming that the linear polarization always 
becomes unstable on account of PS first. What is more, we see that PS is closer to HT 
for low values of s x . We thus expect this situation to be the most favorable to achieve 
squeezing via optical bistability. We checked that varying the parameters Sq and A does 
not change the conclusion. 

To conclude this section, we would like to point out that bistability, as well as PS, 
may disappear when the saturation is too high, as can be seen from FigEJ However, 
we will focus on the low saturation case in the large detuning limit which is the most 
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favorable case for squeezing, and provides analytical results, as well as a clear physical 
understanding. 

4. Mean field fluctuations 

Since we are interested in the quantum fluctuations, we linearize the quantum operators 
around their steady state values following the standard linear input-output method 
[21] . The elliptically polarized solutions are not of great interest for squeezing since the 
predominant circular component sees no atom and the other has negligible intensity. 
Therefore, in all the following, we focus on the linearly polarized state and study how 
both the mean field A x and the orthogonal vacuum field A y may be squeezed. We have 
calculated the outgoing fields noise spectra via a full quantum treatment (see e.g. [23] ) 
involving the four-level system. The outgoing fields are standardly defined from the 
input-output relation [24*j : 

A™l = VTA x , y -A% y (29) 

Yet, to provide clear interpretations as well as analytical results, we derive simplified 
equations, first for the mean field mode A x , then for the vacuum mode A y . 
A similar equation to ()23[) can be derived for the field A x with a term arising from SR 
in SsrA v . In the linearization, this product of zero mean value operators vanishes, so 
that we only have to take saturation into account to derive the spectra of A x . Field 
squeezing owing to optical bistability has been widely studied [213 EE I2Z] and is known 
to occur on a frequency range given by the cavity bandwidth k. The most favorable 
configuration is the bad cavity limit: n is greater than 7 (in our experiment, k ~ 27). 
In the large detuning limit, A ^> g\A x \ ^> 7, the equation for A x reads at order 3 in 
7/A, 

Id s x 
--j-SA x = -(1 + i8 c - iS )SA x + i5 ^- 
ndt 2 

where A x is short for (A x ). This simplified equation yields the classical Kerr terms 
in A 2 X 5A X producing squeezing. Note that absorption, dispersion and the associated 
atomic noise are not included in ([3*U]) . The spectra taking absorption and dispersion 
into account can be easily derived and are shown on Fig[7] The associated susceptibility 
and correlation matrices, [x]xerr and [ajxerr, of the linear input-output theory are 
reproduced in Appendix, and the comparison with a Kerr medium is discussed in [2*5] . 
The situation is more complex for the orthogonal mode on account of SR. 

5. Vacuum fluctuations 

As mentioned in Sec 13.31 SR seems to be a very promising candidate for generating 
vacuum squeezing. However, a careful analysis of the atomic noise, which cannot be 
neglected, is necessary in the squeezing calculations. In the optical pumping regime the 



25 A x + r%SAi 



+ 



A™ (30) 
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Figure 6. Minimal (left) and maximal (right) outgoing field spectra obtained by 
considering only SR effects. Parameters: So — 5, A = 4O7, 7/7.1 = 3, K — 2j, 
s x = O.f , S c = 3. 



circular birefringence, Ssr = —25 J Z /N, is proportional to the ground state population 
difference J z = (022 — 011 )/2 (see Sec IH.4j) . The SR effect is thus closely related to the 
fluctuations of J z , and consequently to the fluctuations of A y via the coupling term in 
A x Ssr [Eq (J2HJ)]- Therefore, we derive general equations for SA y and SJ Z in the Fourier 
domain, and examine their low and high frequency limits. For the sake of simplicity, 
absorption and linear dispersion, again, are not shown; however, the additional terms 
are included in the Appendix. As previously we place ourselves in the large detuning 
limit with s x <C 1 and obtain, discarding terms of order greater than (7/A) 3 , 



— i(u I K)SA y = — (1 + i5 c — i5o)5A y — i5 



A 2 

35 A v - t-^SAI 



+ P(u;)'^A x SJ z + ^=SAf + F Av 



-iu8J~ 



5J Z — 



~ ~2 



N SS Z 

2 Afr, 



+ F Z 



where S z = i{A x A\ — A*A y ) is the usual Stokes parameter (see SecEl) and 



a (a?) 



1UJ 



4 7j 



A(«) 



2 7 — iuj 



ILO 



(31) 
(32) 

(33) 



\{u) ' ^ v ' 4A(w) ' v ' 2( 7 

with Fa v and F z the Langevin operators associated to A y and J z after the 
adiabatical eliminations, 7 P = 7j_s x the optical pumping rate. In (|«U|) . the last term 
of the first line is a crossed Kerr term and clearly contributes to squeezing. We also 
see that the coupling with J z is strongly frequency-dependant and requires a careful 
investigation. In the next sections, we discuss the low and high frequency limits to 
further simplify the previous equations and give simple interpretations for the squeezing. 



5.1. Low frequency: SR effect 

Again, to stress the effects on the fluctuations only due to SR, we neglect the terms 
in s x responsible for the Kerr effect, which will be studied in the next section, and 
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place ourselves in the optical pumping regime, keeping only terms of order 1 in 7/ A. 
Note that this approximation consists in adiabatically eliminating the optical dipoles 
and neglecting the excited state populations and thus limits the analysis to the range of 
frequencies uj <C 7, 7j_. Under these conditions, one has a{yj) ~ 1, /3(u>) ~ 1, \(uj) ~ 1 
and Eq (jH2]) reduces to the linearized optical pumping equation 



— iuoSJ z = — 7 P 



SJ Z - 



N SS Z 
2 I A x \ 2 



+ F Z 



(34) 



It is clear that the fluctuations of J z are governed by the time constant j p , 
consistently with the optical pumping approximation 7 p <C 7. SR is effective only 
at low frequency. Plugging (J54*j) back into (|3*T|). one gets 



^1 + i5 c — iuj/n)5A y = i5 



1 - 



+ i5 Q 



IUJ 

Al 



5A, 



SAT + Fa 



(35) 



'j p - iuj \A X \ 2 

The SR term comes with an amplitude Sq (oc 7/ A) around zero frequency, which is 
much greater than the usual third order saturation non-linearity. Very good squeezing 
could be expected if it were not for the noise coming from the atoms Fa, which we now 
study. The fluctuation operator arising from atomic and field fluctuations reads 

f, 25 A X F z 



N 



+ F A 



(36) 



Jp-lUJ 

The second term Fa is responsible for the noise due to absorption. The first term 
includes the optical pumping noise. One calculates the correlation function of F z via 
the quantum regression theorem ^B! 



(F z (uj)F z (lu')) = 2n5(iu + u')N lp /2 



(37) 



so that 



(F Ay (u) F \M)) 



N 2 rf+LU 2 

,xC 7 



{F z {u)F z {u')) 



2tt5(u 



UJ 



4 7± 7p + 



LL- 



(38) 
(39) 



in which we introduced C = g 2 N/T r y the cooperativity parameter quantifying the 
strength of the atom- field coupling via the cavity (C ~ 100 in our Cs experiment). 
For uj 7 P the noise is thus much more important than the losses due to absorption 
and therefore has a dramatic influence on the squeezing that could have been produced 
by the SR term. Following the method given in [23], we derive the susceptibility and 
correlation matrices which are given in Appendix. We can then calculate the outgoing 
vacuum field spectrum for all the quadratures. Minimal and maximal spectra are plotted 
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on Fig [H] in the "close-to-bad" cavity limit (k = 27) corresponding to our experimental 
configuration. Whereas the first is close to the shot-noise level, the second is extremely 
noisy In the good cavity limit, the noise is even more important. The conclusion is that 
the optical pumping process adds too much noise at zero frequency for SR to generate 
vacuum squeezing. However, this low frequency noise does not prevent squeezing at 
higher frequencies. 



5.2. High frequency limit: crossed Kerr effect 

If one repeats the previous calculation keeping the first order saturation terms in s x and 
considers frequencies uj 3> 7, one finds 

Ns x 6S z (u) , . 

SJz{uj) = T7 W ( } 

This is not surprising, since the evolution times considered are small with respect 
to the atomic relaxation time. The system behaves as if o~ + and <r_ were independent. 
In fact, let us consider two independent two- level systems, 1-4 and 2-3, each with N/2 
atoms. In the large detuning limit, one has 044 = aus + /2 and 033 = a 2 2S-/2, and the 
atomic fluctuations follow the field fluctuations [23] (still at order 3 in (7/ A)) 

r / \ r / \ N 5s + N5s_ 
da u (u) = -da M (u) = -y— , da 2 2{UJ) = -da 33 (u) = -y— (41) 

so that, using 1^4^ | 2 <5(s + — s_) = s x SS z , we retrieve PU|). This equation shows that 
the fluctuations of J z are only caused by saturation and their contribution adds to the 
crossed Kerr terms already mentioned in (|31|) to retrieve a similar "Kerr" equation for 
A y to that of A x at high frequency 



i{uj/K)8A y = -(l + i5 c - i5 )5A 



■ r- S X 

^S - 



2 



25 A v ~ TJtM 



+ + F Ay (42) 



This high frequency behavior is thus characterized by the same Kerr-induced 
optimal squeezing on both polarization modes, consistently with the previous analysis 
for two independent two-level systems. More precisely, the optimal squeezing spectra 
are the same for each mode, but involve orthogonal quadratures [because of the sign 
difference in the Kerr terms between (|3Up and ()42|) ]. We now plot the outgoing fields 
A° ut y squeezing spectra and discuss the squeezing optimization. 

5.3. Squeezing spectra 

To derive spectra for the whole frequency range, we combine both effects by adding the 
matrices obtained in the two asymptotical regimes studied previously. We write the 
complete susceptibility and correlation matrix under the form = [xi^^Kerr + 

[x(c<;)]sfl and [cr(a>)]y = [cr(uj)]Kerr+ [&(u))sr, where the Kerr matrices are those obtained 
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Figure 7. (a) Analytical minimal spectra for A° ut given by the SR effect (dashed), by 
the Kerr effect (light) and by both effects (dark), (b) Exact spectra for the mean field 
mode A° ut (light) and the orthogonal vacuum mode A° Mt (dark). Parameters values: 
S Q = 5, A = 20 7 , 7 / 7 j_ = 3, k = 27, s MI - 0.1, <5 C = 4.6. 

in the high frequency limit, and the SR matrices those obtained at low frequency (see 
Appendix for analytical expressions). This approximation is good since Kerr effect is 
negligible at low frequency compared to SR, while SR breaks down at high frequency. 
In Fig El typical spectra for a working point close to the PS threshold are represented. 
The parameters are chosen to be as close to the experimental situation as possible 
[TT] . We compared these approximate spectra (a) with those obtained with a full 4- 
level calculations based on the linear input-output theory (b). The analytical spectra 
combining Kerr and SR effects show indeed an excellent agreement with the exact 
calculations, as long as the saturation is low. As shown previously, the SR spectrum is 
close to the shot-noise level. The Kerr spectrum is accurate for A y only for u ^> 7 and 
extends on a range of several k as expected |2HJ - The combined spectrum (a) confirms 
that SR destroys completely the squeezing at low frequency and reproduces well the 
exact behavior (b). The best squeezing for A y , obtained at intermediate frequencies, is 
about 25%. 

Note that the Kerr spectrum in (a) is also valid for field A x for all frequencies and 
45% of squeezing is obtained at zero frequency. The situation for the mean field A x is 
identical to that of a circularly polarized field with intensity s x /2 interacting with N/2 
two-level atoms as in for which the Kerr spectrum shows good agreement with the 
exact spectrum. 
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Figure 8. Minimal spectra for the vacuum field (plain) and the mean field (dashed) 
for different working points. 5q — 5, A = 2O7, k — 2j — 67 j_, s max =0.1. The working 
points coordinates are: A (<S C = 5 SR = 4.6, s x /2 = 0.065), B (S c = 5.78, s x /2 = 0.04), 
C (6 C = 6.79, s x /2 — 0.02). The inset shows the working points positions on the 
resonance curve. 



As mentioned in Sec ESI we expect squeezing to improve in the vicinity of the PS 
threshold. We verified this behavior by plotting on Fig |H1 the evolution of the spectra 
when the cavity is scanned while keeping the incident intensity (ex s max ) constant. It 
appears clearly that the best squeezing is obtained at the peak of the resonance curve, 
right before the switching. This is due to the fact that, in the low saturation regime, 
the PS threshold is close to the point where saturation process is the most efficient. 

We then study the effect of saturation and plot on Fig |§] various spectra 
corresponding to working points close to PS with increasing saturation. The conclusion 
is that, for given values of the detuning A and linear dephasing Sq, there is an optimal 
value of s x for squeezing. This is due to the fact that the range for which SR adds 
noise increases with the saturation and eventually destroys Kerr-induced squeezing. 
The optimal saturation value thus corresponds to a compromise between added noise 
and Kerr squeezing in the intermediate frequency range. Therefore, a bad cavity is 
preferable (k ^> 7), since Kerr-induced squeezing occurs on a frequency range given 
by k and SR destroys squeezing for frequencies smaller than 7. Spectra for different 
cavities are represented in Fig E3 The case k = 27 corresponds to the experimental 
situation, " close-to-bad cavity" , the other curves to increasingly bad cavities. Since SR 
is effective on a range smaller and smaller compared to the cavity bandwidth, its effect 
becomes negligible, and 75% of squeezing can be obtained. 

The conclusion is that very interesting squeezing values can be reached in the bad cavity 
limit for both the mean field mode and the orthogonal field mode. In the next section, 
we establish the link between polarization squeezing and the vacuum squeezing obtained 
in our system. 
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Figure 9. Minimal spectra for the vacuum field for different saturation values s max . 
For each value of s max , the working point is chosen close to PS. Parameters: Sq = 5, 
A = 2O7, k = 27 = 67X. 




Figure 10. Minimal spectra for the vacuum for different p = fc/7 (short dash: p = 2, 
long dash: p = 10, plain: p = 50). For each value of /?, the saturation is optimized. 
Parameters: Sq — 5, A — 2O7, 7/71 = 3. 



6. Polarization squeezing 

The noise of the mode with orthogonal polarization with respect to the mean field is 
commonly referred to as polarization noise. However, the study of the polarization state 
fluctuations requires the introduction of the quantum Stokes operators [TJl [TBI EEj 

S (t) = A\A X + A\A y , S x (t) = A\A X - A\A y (43) 
S y (t) = A\A y + A\A X , S z (t) = i{A\A x - A\A y ) (44) 

To be consistent with the definition of our slowly- varying envelope operators A x , 
A y , these Stokes operators are time-dependent and expressed in number of photons per 
second j!4j . They obey the following commutation relationships 



[S (t), Si(t')] = and [Si(t), S 3 (t')] = 2ie ijk S k 5(t - 1') (45) 
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with i,j, k = x,y,z and then the spectral noise densities of these operators, defined 
by (Si(u)Si(uj')) = 2tt5(u + lu^Vs^uj), satisfy uncertainty relations 

VsMVsM > e l:jk \(S k )\ 2 (46) 

The coherent polarization state correspond to the case where both modes A x and A y 
are coherent states. Then the noise densities of the Stokes parameters are constant and 
all equal to Vs^uj) = (So) = \ (A x ) \ 2 + | (A y ) \ 2 for i = 0, x, y, z. The so called polarization 
squeezing is achieved if one or more of these quantities (except Vs ) is reduced below 
the coherent state value 



So 



<1 (i = x,y,z) (47) 



If the mean field is x-polarized, then (S ) = (S x ) = \(A X )\ 2 and (S y ) = (S z ) = 0. 



At first order in noise fluctuations, 8S y and SS Z read 



5S y = \(A X )\ (<L4je*« + 5A y e~^) = \(A X )\ 5X Ay (9 x ) (48) 
5S Z = i\(A x )\ (5Ale ie * - bA y e~^) = \(A X )\ 8X Ay (6 x + tt/2) (49) 

where 6 X is the phase of the mean field and Xa v {0) = A y e td + A y e~ td is the 
quadrature with angle 9 of the orthogonal mode. Therefore the fluctuations of these two 
Stokes parameters are proportional to the quadrature noise of A y and the polarization 
squeezing of S y and S z is simply related to the vacuum squeezing that we have studied 
in the previous sections. The physical meaning of this result is clear: let us choose 
9 X = 0, then geometric jitter on the polarization is due to the intensity fluctuations of A y 
(oc 8S y ), whereas the fluctuations of the ellipticity are caused by the phase fluctuations 
(oc 5S Z ). In the general case, the squeezed and antisqueezed Stokes parameters are 
found to be 



S sq = cos(9 x - 9 sq )S y + sm{9 x - 9 sq )S z (50) 
Santisq = sin(9 x - 9 sq )S y - cos(9 x - 9 sq )S z (51) 

Note that, unlike El HI EI, there is no need to lock the phase-shift difference 
9 X — 9 sq , since it is automatically done in this system; this property appears clearly in 
Eq (|4*2j) where (A 2 X ) /\(A X )\ 2 = e 2ldx . Since the new set of the Stokes parameters So, 
S x , S sq and S antisq still satisfy the relationships (pftjj). we obtain polarization squeezing 
in our system as soon as any quadrature of the vacuum field A y is squeezed, and the 
results of the previous sections can be applied to the squeezed Stokes component. 

7. Conclusion 

We have presented a study of polarization switching in an X-like 4-level atoms ensemble 
illuminated by a linearly polarized light in an optical cavity. PS has been traced to 
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self-rotation and simple criteria allow for a clear understanding of the switching effects 
and the multistable behavior of the system. The steady state analysis enables one to 
figure out the interesting working points for squeezing. 

In terms of squeezing the respective contributions of SR and saturation have been 
investigated and compared to a full quantum calculation. Since the propensity for 
squeezing of SR is cancelled by atomic noise at low frequency, the squeezing originates 
from Kerr effect. The mean field mode is squeezed via the usual saturation effects, 
whereas the vacuum mode squeezing is induced by the mean field via crossed Kerr effect. 
Both SR and crossed Kerr effects can be dissociated in a bad cavity configuration, thus 
allowing for high squeezing values. Last, this vacuum squeezing is shown to be equivalent 
to squeezing one Stokes operator. 

Appendix A. 

Using the input-output theory notations we give here the expressions of the 

susceptibility matrix [x(<^0] and the correlation matrix [cr(u;)] for field A y . In the high 
frequency limit, they resume to those derived in [23] in the large detuning limit 



e — — 1 yields the susceptibility matrix for the vacuum mode. To retrieve the matrix 
for A x , e should be taken equal to +1. This matrix corresponds to approximating the 
atoms ensemble with a Kerr medium: the term of order 1 in 1/ A is the linear dephasing, 
the second order matrix represents dispersion and absorption and the third order term 
is the non-linear dephasing corresponding to the Kerr effect. The associated correlation 
matrix is 



In the Kerr limit, the atomic noise comes only from the frequency independent 
linear losses of the Kerr medium, which acts as a beamsplitter for the field. Similar 
matrices can be derived for field A x in agreement with (jHOjl . At low frequency, however, 
the previous matrices have to be completed by 




(A.l) 




(A.2) 




(A.3) 
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I ^ f ^ + ^ P \ rA4 x 

For u; 7 p , the vacuum correlation matrix is equivalent to 

so that a lot of noise is reported on all the quadratures of A y for frequencies of the 
order of j p , as pointed out in Sec 15.11 For frequencies u 3> 7, the SR noise terms vanish, 
allowing for crossed Kerr effect to produce squeezing. 
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